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Abstract 

P. W. Anderson raised an important question in 2007: Is There Glue in Cuprate 
Superconductors? The author believes that the change of the electron clouds of ions is the 
glue in cuprate superconductors. The change of the electron clouds of the ions in the parent 
structure of the layered high-temperature superconductors CaCuO2 has been studied by the 
first-principles calculations. The electron clouds of Cu” and O” ions change obviously under 
electric fields. It is also found, for the first time, the characteristic frequencies of the change 
of the electron clouds are 250 meV, 360 meV, and 100 meV, respectively, for the modes 
observed. The frequencies are low and close to that of lattice vibrations, indicating the change 
of the electron cloud of ions can be the electron-pairing medium in cuprate superconductors. 
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1 Introduction 

In 1986, J. G. Bednorz and K. A. Müller !!! discovered that there may be 
high-temperature superconductivity (HTS) in La-Ba-Cu-O oxides. After that, many cuprate 
superconductors were discovered with transition temperatures higher than 77K PPI, The 
electron-lattice interaction cannot explain the electron pairing mechanism in cuprate 
superconductors. The pairing mechanism for unconventional HTS is still under debate. 

P. W. Anderson "! raised an important question in 2007: Is There Glue in Cuprate 
Superconductors? The author believes that the change of the electron clouds of ion can be the 
glue in cuprate superconductors. The mechanism is as follows. When a free electron comes to 


a new location, the electron cloud of the ions in the vicinity will change. In this way, the 
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charge density around the free electron will decrease. When the free electron leaves, the 
electron cloud of the ions will not relax immediately, so that there will be a region lack of 
charge (equivalent to a positive charge region) that will attract another free electron. 
Attraction between free electrons appears. This mechanism is essentially the same as the 
electron-phonon interaction, except that the medium is the change of the electron cloud, not 
the displacement of ions. 

According to Born Oppenheimer approximation '°!, the electrons are moving and 
responding to forces very quickly, because the electrons have much smaller masses than the 
nuclei (more than 1000 times). The electron density changes too fast and can not be excited 
by free electrons. So, it is generally believed that electron pairing cannot be achieved by the 
change of the electron cloud. 

But, can the electron density change as slowly as the nucleus? If it is, then it can be 
the electron-pairing medium. 

Based on the above consideration the author made an investigation. The research object 
is the parent structure of the layered high-temperature superconductors CaCuO, 1, Firstly, the 
change of the charge density under the influence of electric fields was studied. Secondly, the 
characteristic frequency of the change was calculated with real-time TDDFT method ® °, The 
author found, for the first time, the characteristic frequencies of the change of the electron 
clouds of ions are close to that of the lattice vibration. So, the change of the electron cloud of 
ions can be excited by free electrons, indicating the change of the electron clouds can be 
electron-pairing medium in cuprate superconductors. This paper will report the methods and 


results. 


2 Methods 

Firstly, VASP!" and QE"! are used to study the change of charge densities under the 
influence of electric fields. In calculations with VASP, the projector augmented plane wave 
(PAW) method I4 was used to describe the interaction between the valence electrons and the 
ion cores. The generalized gradient approximation (GGA) of Perdew-Burke-Ernzerhof (PBE) l5! 
was used for the exchange correlation potential. To describe the strong correlation of the 3d 


electrons, the GGA+U method is used. For the 3d orbitals of Cu, the Uef value is 6.5eV HIN, 
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Plane waves were used to expand the valence electron wave functions. The cutoff energy is set to 
400 eV which has been tested by taking both accuracy and computing cost into account. The 
valence electron configurations of the pseudopotentials are taken as 3s73p°4s for Ca, 3p°3d'°4s! 
for Cu and 2s?2p* for O. Further calculations with the Heyd, Scuseria, and Ernzerhof (HSE) "®! 
hybrid-functional were performed for cross validation. In the calculation with QE, 
norm-conserving pseudopotentials were employed to calculate the change under uniform electric 
fields. 


Then, the characteristic frequency of the change of the charge density was studied by the 


real-time TDDFT method. Calculations were conducted within the Octopus package l” *"!, 


3p°4s* of Ca, 3p°d'°4s' of Cu and 2s’p* of O are taken as valence electrons. The generalized 


gradient approximation (GGA) of Perdew-Burke-Ernzerhof (PBE) was used to describe the 
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exchange-correlation energy. The GGA+U method | was considered to deal with the 
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strong correlation of the 3d electrons. HSCV pseudopotentials ! were adopted. 


Approximated Enforced Time-Reversal Symmetry (AETRS) algorithm was used to 


approximate the evolution operator and the time step is 0.002 ii/eV. 


3 Results and Discussions 

Fig. 1 displays the crystal structure and the charge density difference (CDD) with and 
without electric fields. The yellow color represents a positive value or increase of the charge 
density, while the blue color negative or decrease. The charge density around Cu** ions 
change significantly. There is no similar change around Ca” ions. The change is more like a 
rigid rotation rather than an elastic deformation, because the change is not entirely along the 
direction of the electric field. Some areas increase, while some areas decrease. Furthermore, 
the pattern of the change is like 3d electron clouds. Why do electron clouds change this way? 
We know that transition metal ions have incompletely filled 3d orbitals. The electron cloud 
does not have perfect spherical symmetry. It is like a polar molecule and will change under 


electric fields. 


The charge density around O” ions changes significantly, too. The author thinks that the 
valence of O* ions is not -2 strictly. If it is a strict -2 valence, it will remain unchanged just 
like Ca? ions. The author has also studied eight other unconventional superconductors 


(Fe2KSe2, LazFe2As2O2, Nd2Fe2As2O2, BazFesAss4, YBarCu307, HgBazCa2Cu3Os, 


TlsBazCaCu20s and Bi2Sr2Ca2Cu3O10) °“. Under static electric fields, the electron clouds of 


transition metal ions change significantly. 


Fig. 1 (a) Crystal structure of CaCuO2; (b) Charge density difference (CDD) caused by Uniform E-field (QE, 
isosurface = 0.005, efield_cart(1) = efield_cart(2) = efield_cart(3) = 0.002, kpoint = 5x5x7, ecutwfc = 60.0 Ry ); (c) 
CDD caused by changing one atom slightly (VASP, GGA+U, kpoints = 6x6x10, isosurface = 0.005); The effect of 

an electric field induced by a slight change in the position of one oxygen atom in the copper oxide plane. The 

original coordinate of the oxygen atom is (0.75, 0.25, 0.00), and changed to (0.8, 0.3, 0.00); (d) CDD caused by 


changing one atom (VASP, HSE, kpoints = 6x6x10, iso = 0.005). Plots were generated using VESTA P5, 


In order to study the real-time evolution of the change of the charge density, two 
methods are used to excite the charge vibration. One is to apply a time dependent electric 
field, while the other is to slightly change the position of some atom. The form of the 


time-dependent electric field is as below: 


t—2t,—- 


FO =F cos? 0) cost If \t—t,|>). then f(t) =0. 


To 


Fig. 2 shows the crystal structure and the real-time evolution of the charge density 
excited by a time dependent electric field with œ = 1.8 eV and Fo = 1.0 eV/angstrom. The 
excitation vanishes after 5000 steps and the charge densities vibrate freely thereafter. After 
7000 steps, the charge densities are recorded. 100 represents the charge density difference 
between the 7100th step and the 7000th step and the same below. The change of the electron 
clouds of the O* ions (marked with an arrow) become obvious gradually with the evolution 


steps. The change of the charge density reaches its maximum after about 500 steps, and the 
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corresponding time is 1.0 f/eV. It is the time from zero to the maximum. The time of one 


period should be 4x1.0 fi/eV and the frequency is about 250 meV. 
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Fig.2 Crystal structure and the charge density evolution of CaCuOz. The isosurface is 0.002 e/bohr?. The yellow 
color represents an increase of the charge density, while the blue a decrease. w = 1.8 eV and Fo = 1.0 eV/angstrom. 


The excitation vanishes after 5000 steps and the charge densities vibrate freely thereafter. 


Fig. 3 shows the crystal structure and the real-time evolution of the charge density 
excited by a time dependent electric field with œ = 1.8 eV and Fo = 10 eV/angstrom. The 
excitation vanishes after 5000 steps and the charge densities vibrate freely thereafter. After 
18000 steps, the charge densities are recorded 50 represents the charge density difference 
between the 18050th step and the18000th step. The change of the charge density of the Cu** 
ion marked with an arrow reaches its maximum after about 350 steps, and the frequency is 


about 360 meV. 
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Fig.3 Crystal structure and the charge density evolution of CaCuO> excited by a time dependent electric field with 
æ = 1.8 eV and Fo = 10 eV/angstrom. The isosurface is 0.04. 

Fig. 4 shows the crystal structure and the real-time evolution of the charge density 
excited by slightly changing an Ca atom. To excite the charge vibration, the position of a Ca 
atom is slightly changed. The original coordinate is (0, 0.5, 0.5), and modified to (0.02, 0.52, 
0.5), ie., the Ca atom was moved in the (001) direction by 0.15 A. Ground state was 
calculated. Then the coordinate was modified back to (0,0.5,0.5) for TDDFT calculation. 
Charge densities were recorded after 13000 steps. The electron cloud of the Cu” ion marked 
with an arrow changes gradually. The change reaches its maximum value after1300 steps and 


the frequency is about 100 meV. 
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Fig. 4 Crystal structure and evolution of the charge density of CaCuOz. The isosurface is 0.03. 


The author also studied the change of the electron cloud of the ions in other 
superconductors °!, For the TDDFT study of La:Fe2As202, FeSe sheet and HgBarCa2Cu3Ox, 
the characteristic frequencies are 160 meV, 190meV, and 250meV, respectively 27 For 
BaFezAs2, the characteristic frequencies are 150 meV, 160 meV, 250 meV, and 200 meV, 
respectively, for the modes the author observed P®!, Table 1 gives the maximum phonon 
frequencies and Te of three typical conventional superconductors. The characteristic 


frequencies of the change of the electron clouds are close to that of the lattice vibrations. 


Table 1 Maximum phonon frequencies and T: of three typical conventional superconductors. 


Pbl] MgB2l°] H3S80l 
w / meV 9 90 250 
Te/K 7 39 164 


It is worth noting that the change of the electron clouds of transition metal ions is very 
complicated. There may be different modes, corresponding to different frequencies. The 
results are unexpected, because the general view is that the change of the electron density is 
very quick and the frequency is much higher than the lattice vibration. The frequencies the 
author obtained are close to that of the lattice vibration, indicating it can be excited by free 
electrons. So, the change of the electron clouds of transition metal ions can be the 


electron-pairing medium. W. A. Little P" 


gave similar results, but the frequencies are too high 
(> 1000 meV). Whether it can be excited by free electrons should be justified. 
Why do the electron clouds of transition metal ions change this way? The main reason is 


that the electron clouds possess no spherical symmetry, and easy to change under electric 


fields. 


4 Summary 
The change of the electron clouds of ions in CaCuO2 was studied by the first-principles 
method. The author found, for the first time, that the electron cloud of ions can change slowly. 


The frequencies of evolution of electron clouds match well with the frequencies of the lattice 


vibrations. This is contrary to Bonn Oppenheimer approximation. Though the frequencies 
obtained by this method are not accurate, it can give a significant evidence that the change of 
electron clouds can be the electron-pairing medium in HTS. The change of electron cloud can 


be used as a pairing medium for superconducting electrons. 
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